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Abstract 

We examine random matrix models of QCD which are capable of supporting both chiral and 
diquark condensation. A numerical study of the spectral densities near zero virtuality shows that 
the introduction of color in the interactions does not alter the one-body results imposed by chiral 
symmetry. A model with three colors has the spectral density predicted for the chiral ensemble 
with a Dyson index f3 = 2; a pseudoreal model with two colors exhibits the spectral density of the 
chiral ensemble with (3 = 1. 
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Chiral random matrix theory (%RMT) is based on the observation that many of the 
low-energy properties of QCD are dominated by its global symmetries [lj]. Random matrix 
models |2| thus attempt to capture the basic mechanisms for chiral condensation by reducing 
the QCD interactions to their essential structure. These models introduce low-lying modes 
which respect a basic left-right symmetry but which interact via random matrix elements. 
The consequences of chiral symmetry have been investigated at two levels. The microscopic 
level deals with the statistical properties of the eigenvalues of the Dirac operator and their 
correlations. Here, %RMT studies 3,0,0] have helped in understanding how these properties 
are determined by the spontaneous breaking of chiral symmetry as indicated, for example, 
by the universality of the spectral density near zero virtuality and associated sum rules [lj]. 
At the macroscopic level, it is possible to consider the consequences of chiral symmetry for 
the global state of the random matrix system. With the aid of the mean field approximation 
and additional prescriptions for including the effects of temperature and chemical potential, 
it is straightforward to construct the partition function for xRMT models and thus examine 



the patterns of chiral symmetry breaking as a function of T, fj,, and the quark mass 
Many of the properties of the resulting phase diagram are direct consequences of chiral 
symmetry and largely independent of the detailed form of the interactions. They are thus 
expected to provide guidance in our understanding of the QCD phase diagram particularly 
in cases (e.g., N c = 3 and /i ^ 0) where the Dirac operator is non-Hermitean. These are 
cases where lattice simulations cannot rely on importance sampling and are very difficult to 
perform [3, S, Z3, EH - 

There has been considerable interest in the possibility that both chiral and diquark con- 
densates can develop and com pet e thermodynamically in QCD. (For a review of diquark 
condensation in QCD, see , ^j.) Thus, we recently proposed a random matrix 



model for QCD which goes beyond %RMT and has interactions which implement coexisting 
chiral and color symmetries We have examined this model at the macroscopic level 

by identifying the allowed topologies for the (T, fi) phase diagram. The phase structure 
is determined by a single parameter a, defined as a ratio between coupling constants in 
the chiral and diquark condensation channels. Physically, this ratio measures the relative 
strengths of the chiral and color symmetries in the random matrix interactions. In the cases 
relevant for QCD [either SU(2) or SU(3)], a has a fixed value and is associated with a given 
phase structure. The topology of this phase structure is robust with respect to moderate 
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variations in the detailed form of the interactions (i.e., variations in a). 

The addition of color structure to the interactions implies additional constraints on the 
random matrix ensemble considered. As we shall see below, the color generators modify 
the statistical weights of the interaction matrix elements. The question then arises whether 
these constraints are capable of altering the various results associated with chiral symmetry 
alone. Our previous considerations suggest that this is not the case at the macroscopic level. 
In the case of an interaction which is completely dominated by chiral symmetry (a —>■ 0), 
our models precisely reproduce the phase diagram of xRMT (Ref. p). In cases relevant for 
QCD with fundamental fermions |2y] and three colors (iV c = 3 and a = 0.75), the diquark 
phase develops in regions where it is thermodynamically advantageous, but color does not 
otherwise modify the phase structure. Thus, for N c = 3, color does not seem to weaken the 
chiral correlations at the macroscopic level. 

For N c = 2, the situation is more subtle: the gauge interaction is pseudoreal, and baryons 
and mesons belong to the same multiplets. For m = 0, \x = 0, and Nf flavors, the quark 
Lagrangian has an extended SU(2iV/) invariance which relates chiral and diquark fields. 
This flavor symmetry is explicitly broken for either m > or // > 0. Studies of the sym- 



metry breaking patterns showed that a 



~ oc m 1 / 2 and T < T c [l9|, |20j, |21|. Our model exhibits the SU(2iV/) symmetry 



diquark condensed phase becomes favorable for 



and produces results which agree with the general predictions of Ref. 18] as well as with 

nn 

those of chiral perturbation theory [TjJ |2Jj. Thus, for N c = 2, the introduction of color in 
the random matrix interactions does not lead to unexpected topologies of the phase diagram 
(i.e. unexpected macroscopic properties). 

The purpose of this paper is to investigate the consequences of the additional color corre- 
lations at the microscopic level. To this end, we examine the distributions of the eigenvalues 
of the Dirac operator near zero virtuality and compare them with the analytic forms pre- 
dicted by xRMT. The rest of the paper is organized as follows: We present the models in 
Sec. HI discuss aspects of their macroscopic and microscopic spectral densities in Sees. ITT1 and 
comment on other properties in Sec. IIV[ and conclude in Sec. El 
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I. THE RANDOM MATRIX MODELS 



To understand the form of the correlations induced by each symmetry, it is useful to 
provide a brief description of the models of Refs. Q], which we will refer to as I. In the present 
paper we will restrict our attention to a theory with Nf light flavors, zero temperature, zero 
chemical potential, and zero quark mass. We start by recalling the basic form of chiral 
random matrix theory, in which only chirality is introduced, and then extend the model to 
include color. 

We consider first a chiral random matrix model and work in the sector of zero topological 
charge for simplicity. The partition function has the form 

ft 



N f 

Z = ( DW J] D^* ex P 

J 1=1 



N 

i 



eX p(-^^Tr[Wt]), (1) 
i=i J A 

where ip i and ip* are independent Grassmann variables representing the quark fields and 
where the matrix T> represents the Dirac operator. Its block structure reflects chiral sym- 
metry Working in a suitable basis of left and right states (1 ± 75)0 n , T> has the form |2[ 



V 



f 

i W 1 j 



(2) 



where W is an N x N block matrix. The integral in Eq. $IJ is over the matrix elements of 
W, DW is a Haar measure, and Tr denotes a trace over the N matrix indices. The model is 
thus a theory of 2N low-lying modes which respect chiral symmetry and whose interaction 
matrix elements Wij are drawn on a Gaussian distribution. The number of modes scales 
with the volume of the system; the thermodynamic limit is taken as N —>■ oo. 

For later comparisons, it is useful to understand how the gauge group is taken into 
account. In a model which implements the global symmetries of QCD with N c = 3, the 
matrix elements of W are complex. For the analysis below, it is worth noting that their real 
and imaginary parts satisfy no particular relationship and are thus drawn independently. 
This case corresponds to a Dyson index (3 = 2 and is described by the chiral unitary ensemble 
(XGUE). For N c = 2, the gauge group „ pseudoreal ™ mentioned above. The Dirac operator 
then contains an additional antiunitary symmetry [2|], which allows one to choose a basis in 
which the matrix elements of W are real. This case leads to the chiral orthogonal ensemble 
(xGOE). In short, the value of iV c enters ^RMT only through the reality of the matrix 
elements of W. 
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In model I, we included color directly in the interactions in a way which mimics single- 
gluon exchange. For two flavors, the partition function takes the form 



DH Dip\ L>Vi Dip* Dift exp 
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T 



T 



( ^ n \ / 




(3) 



where ip\ and ipj denote the quark fields for flavor 1 and flavor 2, respectively, DH is a 
measure to be defined below, and T> c is the single-quark Dirac operator. Note that the 
subblock associated with flavor 2 has been transposed in order to exhibit the possibility of 
forming (ipjiftx ) condensates (see I). 

The diagonal block T> c now reflects both chiral and color symmetries. In order to be 
able to define an order parameter which is antisymmetric under the permutation of an odd 
number of quantum numbers, spin has to be introduced together with color 16]. The Dirac 
operator T> c then has the chiral structure of Eq. (j2j) where W is exploded into a 2N C x 2N C 
matrix of embedded spin and color subblocks, 

3 A^-l 

w = H E K) V- (4) 

/i=0 a=l 

Here cr M = (l,ia) are 2x2 spin matrices, A a are the N c x N c matrices of SU(iV c ) and A^ a 
are n x n real matrices representing the gluon fields. Taking into account all substructures, 
each matrix W is thus N x N with iV = 2N c n. The matrix elements of A^a in Eq. (j3J are 
distributed according to the measure DH, Eq. (jSJ), which takes the form 

DH = III DAA exp (-NT? J2 Tr[A M a (A-fl) > (5) 

where DA^ a are Haar measures, So is a constant, and the superscript T denotes a transpo- 
sition. Again, the thermodynamic limit corresponds to iV — ► oo. 

Evidently, these spin and color substructures impose strong constraints between the real 
and imaginary parts of the matrix elements of W. Consider, for example, the contributions 
to V c , Eqs. (121) and of the random matrices A i (/x = and a = 1) and A 02 . Their 
matrix elements are real and are drawn independently. When combined to form W as 
prescribed in Eq. (0J), the matrix elements of Aoi multiply ero, a real diagonal spin matrix, 
and Ai, a real color matrix. Hence, Aoi contributes to the real part of W only. Similarly, the 
matrix elements of A 02 combine with <jq (real) and A2, an imaginary color matrix. They thus 
contribute to the imaginary part of W . Hence, in contrast to ordinary xRMT, the matrix 



elements of W are complex for both N c = 3 and N c = 2. The real and imaginary parts 
of W arise from well-defined combinations of the matrix elements of A^ a . Their statistical 
distributions are then dictated by the content of the spin and color block matrices, which 
thus introduce well-defined correlations. 

Do these additional correlations preserve those imposed in Eq. (J2J) by chiral symmetry? 
The W matrix elements are complex for all N c ; their real and imaginary parts are no longer 
independent random variables. One might thus be concerned that the statistical properties 
of the eigenvalues of the Dirac operator would differ from those of the xGUE and %GOE. 
We now consider a number of spectral properties to indicate that this is not the case and 
that the additional color symmetries do not alter the statistical properties due to chiral 
symmetry. 

II. THE MACROSCOPIC SPECTRAL DENSITY 

As an initial measure of the statistical properties, we consider the spectral density p(A), 
defined as 



where Aj are the 27V" eigenvalues of the Dirac operator and () denotes an ensemble average. 
While p(A) is not universal in the usual random matrix sense, we will focus on those symme- 
try properties which are expected to be protected. Consider first the chiral random matrix 
model defined in Eq. (JTJ). Because of the block structure of V, Eq. (j2J), the eigenvalues Aj 
occur in pairs of opposite signs. Moreover, Eq. (0) shows that the spectrum is iVy-fold de- 
generate. Consider next model I, Eq. Q, which may seem different at first glance because of 
its more elaborate structure. In fact, each flavor subblock has the same chiral substructure 
as in Eq. (J2J), and the eigenvalue spectrum remains symmetric about A = 0. Further, the 
eigenvalues of — T>^, Eq. (j3J), are degenerate with those of V c so that the spectrum is again 
twofold degenerate (Nf = 2). Hence, the same basic chiral and flavor symmetries prevail in 
the two models. 

In order to facilitate numerical evaluation, we now consider the quenched limit Nf — > 0. 
This limit is free from contributions from vacuum graphs (through powers of the determi- 
nant of the Dirac operator in the partition function; see ]]]). The ensemble average of a 




(6) 



6 



given quantity then amounts to a mere counting of the contributions from the individual 
eigenvalues of V (or V c ), distributed according to the normal laws of Eq. (JTJ) or Eq. ©, as 
appropriate. 

For large matrices, we find that the model I is numerically consistent with the (non 
universal) semicircle law familiar from the xGUE and xGOE, 



, .. 1 - (^) if |A| < 2/E, 
Jim p(A) = M V 2 J M - I (7) 

otherwise. 

In the case of chiral random matrix models, E is the variance of the distribution in Eq. ((TJ) . 
In model I, £ is proportional to E , Eq. (|5|). in a manner that will be discussed shortly. Note 
that thanks to the Banks-Casher relationship [22], £ is in all cases to be identified with the 
chiral order parameter: 

(ipip) = lim lim np(\) = S. (8) 

Having obtained the semicircle distribution in model I is a natural result. The random 
matrix interactions of Eq. (J3J) mix a set of 4 x (Af? — 1) independent real matrices, A^, in 
a democratic way. This ensemble naturally leads to the semicircle law familiar from most 
elementary random ensembles, including the chiral ones. 

Further remarks can be made about the dependence of E with respect to N c . For xRMT, 
E does not depend on (3, and hence does not depend on N c . In model I, however, E is a 
function of N c and Eo which can be easily determined by noting the following relationship 
between the radius of the semicircle 2/E and the variance (TrjWiy"!]): 

(TrpWt]> = ^. (9) 

From the definition of W, Eq. (j3J), and the distribution in Eq. (j5J, we have 

o(N 2 — 1 1 N 

(Tr[Wt]) = ±Y,(n\ a (A, a f}) = [ C N2 } - 2 , (10) 

which gives 



P {N 2 C - 1) 



0.75 En if N R = 3 



0.8165 E if A^ c = 2. 



(11) 
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The N c dependence of the chiral condensate reflects the fact that the strength of the inter- 
actions in the chiral channel varies with N c . In fact, the ^-dependent prefactor in Eq. (fTU|) 
is directly proportional to the Fierz coefficient which appears when the random matrix 
interactions are projected onto the chiral condensation channel (Ref. Q]). 



III. THE MICROSCOPIC SPECTRAL DENSITY 



The structures that are characteristic of chiral symmetry are most clearly seen in the dis- 
tributions of the eigenvalues near zero virtuality. It is thus worth considering the microscopic 
spectral density, defined at the level of single eigenvalues as 

p s {z) = lira — p , (12) 



n-ico 2iV£ r V2AE. 

where £ is the chiral order parameter of the theory at hand. The microscopic spectral 
densities of the various chiral ensembles are strikingly different. In the quenched limit, 
xGUE gives 



while xGOE gives 



23 



24] 



Ps(z) = - (J 2 (z) + Jl(z)) ( X GUE, N f = 0, z > 0), (13) 

Q 

Ps{z) = \ (Jliz) - Hz) J 2 (z)) - i J (z) Qf*d« J 2 (u) 

( X GOE,N f = 0,z>0). (14) 

Figures 1 and 2 show the histograms obtained in the models of Refs. [16], plotted against 
Eq. ffTTTf) for the model with N c = 3 and against Eq. (fPfj) for that with N c = 2. The data 
were obtained from an ensemble of 2 x 10 5 matrices with blocks W of size 120 x 120 (i.e., 
2nN c = 120). As in xRMT, the behavior of the small eigenvalues depends strongly on 
whether N c = 3 or N c = 2. In fact, the data agree with the theoretical curves for the 
corresponding chiral ensembles to expected statistical accuracy. The microscopic spectral 
density thus appears to preserve the form dictated by chiral symmetry in spite of the presence 
of additional color symmetries. In particular, we observe the familiar depletion near z = 
for N c = 3. 

The spectrum in Fig. 2 is flatter than that for N c = 3. This is a consequence of the 
greater motion of the individual eigenvalues in the model with N c = 2. This difference is 
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familiar in chiral ensembles, where it is interpreted as a consequence of the fact that the 
spectrum of a real matrix (xGOE) is less rigid than that of a complex matrix (xGUE). It is 
interesting to see that the level motion in our model with N c = 2 agrees with that expected 
for xGOE even though the matrix elements of W are no longer real. Again, the correlations 
induced by chiral symmetry are maintained in spite of the additional constraints introduced 
by the color structure of W. 

Of course, the preservation of the N c = 2 microscopic spectral density in spite of the 
complexity of W should not come as a surprise. The interactions in Eq. Q are, in fact, 
pseudoreal when N c = 2 [l^. Thus, we could have used an appropriate basis of states 
such that all matrix elements have vanishing imaginary parts. The model would then have 
displayed a greater similarity to a xRMT model with (3 = 1, and it would thus be natural 
to find a spectrum which reproduces that of xGOE. 



IV. OTHER SPECTRAL PROPERTIES 



We have also studied short range eigenvalue correlations by examining the distribution 
of level spacings, p(s). Here, s is the nearest neighbor spacing measured in units of the 
local average spacing. (Hence, the average value of s is 1 by definition.) Note that p(s) is a 
bulk observable which is sensitive to chiral symmetry only for small s. Figure 3 shows the 
relative difference between the results obtained for the model of I with N c = 3 and those of 
the xGUE. The data were obtained from a set of nine independent series of diagonalizations 
of 10000 matrices of size N x N = 120 x 120. To avoid side effects due to the finite 
range of the random matrix support, we only kept a central portion of the spectrum (30th 
eigenvalue to 90th). The level spacing distribution for a given model is taken as a bin- 
by-bin average over the nine series of runs, with error bars corresponding to the variance 
in each bin. The plot shows Sp(s) = [p x gve(s) - Pn c =3(s)]/p xGVE (s), where £> x gue(s) is 
the level spacing distribution for xGUE and Pn c =3{s) is that for model I with N c = 3. 
[Note that p x gue(s) = Pgue(s).] The level spacing difference is everywhere consistent with 
zero. Similar results were obtained for a comparison with xGOE. Thus, the level spa cing 
distribution reveals no statistically significant differences between the models of Refs. (lfij 
and their corresponding chiral ensembles. 
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It is also useful to construct the variance-covariance matrix D, defined as 



A J = ((A,-A l )(A,-A,)), (15) 

where \ is the mean value of eigenvalue i. Numerical studies of the models of Refs. 
indicate that all eigenvalues of D are greater than zero. Thus, the eigenvalues of the Dirac 
operator do not satisfy any linear constraints. The eigenvalues of D are related to the sta- 
tistically independent fluctuations of the eigenvalues of the Dirac operator about their mean 
values. (This relation would be exact if the joint probability distribution for the eigenvalues 
of the Dirac operator were strictly Gaussian.) Numerical studies reveal good agreement 
with the spectrum of normal modes known analytically for the xGOE and xGUE [25]. 
This provides an additional indication that the inclusion of color symmetry has not altered 
the statistical properties of the models of Refs. from those of the corresponding chiral 
ensembles. 



V. CONCLUSIONS 



We have investigated a number of properties of the eigenvalue spectrum of random ma- 
trix models which include both chiral and color symmetries. We find no deviation from the 
analytic results of %RMT for either the microscopic spectral density or the level spacing 
distribution. Given the relatively elaborate block structure of these models, such studies 
are most easily performed numerically. While more complicated spectral correlators have 
not been investigated, we find no grounds to doubt that they will also reflect the underlying 
chiral structure of the problem. (It should be emphasized that this extended block structure 
does not lead to complicated forms for the partition function and that studies of the macro- 
scopic properties of these models are straightforward for all T and fi.) The spin and color 
block structures of the interactions do not seem to upset spectral features associated with 
chiral symmetry. Rather, these additional correlations appear to act in a channel which is 
"orthogonal" to the chiral channel. We do not regard this result as surprising. One may 
think of the models of Refs. as schematic lattice calculations which implement single- 
gluon exchange between discrete quark fields. This "lattice calculation" does not possess 
any other length scale than that introduced by £ and is, hence, free of a Thouless energy. 
Because chiral symmetry is respected by the interactions, the results of %RMT should be 
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anticipated at all small energies. This is indeed what our numerical analysis reveals. 

The inclusion of color does, of course, introduce additional correlations in the spectrum 
of the Dirac operator. It would be of interest to construct new spectral measures which 
could probe these correlations. In particular, it would be valuable to find a color analogue 
of the Banks-Casher relation to facilitate the determination of the strength of the diquark 
condensate. 

Acknowledgments 

We thank J. J. M. Verbaarschot for provocative questions which initiated this work and 
D. Toublan and K. Splittorff for useful discussions. 



[1] For a review of xRMT, see J. J. M. Verbaarschot and T. Wettig, Annu. Rev. Nucl. Part. Sci. 
50, 343 (2000). 

[2] E. V. Shuryak and J. J. M. Verbaarschot, Nucl. Phys. A560, 306 (1993); J. J. M. Verbaarschot, 
Phys. Rev. Lett. 72, 2531 (1994); Phys. Lett. B 329, 351 (1994). 

[3] J. J. M. Verbaarschot and I. Zahed, Phys. Rev. Lett. 70, 3852 (1993). 

[4] J. J. M. Verbaarschot, Nucl. Phys. B427, 534 (1994). 

[5] J. J. M. Verbaarschot, Nucl. Phys. B426[FS], 559 (1994). 

[6] A. D. Jackson and J. J. M. Verbaarschot, Phys. Rev. D 53, 7223 (1996). 

[7] M. A. Halasz, A. D. Jackson, R. E. Schrock, M. A. Stephanov, and J. J. M. Verbaarschot, 
Phys. Rev. D 58, 096007 (1998). 

[8] M. G. Alford, A. Kapustin, and F. Wilczek, Phys. Rev. D 59, 054502 (1999). 

[9] S. Chandrasekharan and U. J. Wiese, Phys. Rev. Lett. 83, 3116 (1999). 
[10] Z. Fodor and S. D. Katz, Phys. Lett. B 534, 87 (2002); J. High Energy Phys. 03, 014 (2002). 
[11] P. de Forcrand and O. Philipsen, Nucl. Phys. B642, 290 (2002). 

[12] B. Barrois, Nucl. Phys. B129, 390 (1977); D. Bailin and A. Love, Phys. Rep. 107, 325 (1984). 
[13] R. Rapp, T. Schafer, E. V. Shuryak, and M. Velkovsky, Phys. Rev. Lett. 81, 53 (1998); Ann. 

Phys. (N.Y.) 280, 35 (2000). 
[14] M. Alford, K. Rajagopal, and F. Wilczek, Phys. Lett. B 422, 247 (1998). 

11 



[15] For a review on color superconductivity, see K. Rajagopal and F. Wilczek, B. L. Ioffe 
Festschrift, At the frontier of particle physics/Handbook of QCD, edited by M. Shifman (World 
Scientific, 2001), Vol. 3, p. 2061 |hep-ph/0011333[ 



[16 

[17; 

[18 
[19 

[20 
[21 
[22 
[23 
[24 
[25 
[26 



B. Vanderheyden and A. D. Jackson, Phys. Rev. D 61, 076004 (2000); 62, 094010 (2000). 

B. Vanderheyden and A. D. Jackson, Phys. Rev. D64, 074016 (2001). 

J. B. Kogut, M. A. Stephanov, and D. Toublan, Phys. Lett. B 464, 183 (1999). 

J. B. Kogut, M. A. Stephanov, D. Toublan, J. J. M. Verbaarschot, and A. Zhitnitsky, Nucl. 

Phys. B582, 477 (2000). 

K. Splittorff, D. T. Son, and M. A. Stephanov, Phys. Rev. D 64, 016003 (2001). 

K. Splittorff, D. Toublan, and J. J. M. Verbaarschot, Nucl. Phys. B620, 290 (2002). 

T. Banks and A. Casher, Nucl. Phys. B169, 103 (1980). 

P. J. Forrester, P. Nagao, and G. Honner, Nucl. Phys. B533, 601 (1999). 

B. Klein and J. J. M. Verbaarschot, Nucl. Phys. B588, 483 (2000). 

A. D. Jackson, C. B. Lang, M. Oswald, and K. Splittorff, Nucl. Phys. B616, 233 (2001). 
In this work, we consider only fundamental fermions. For models with adjoint fermions, see Q, 



19 



20, 



21] 



[27] This can be established by considering the resolvent operator G(z) = (Tr(z — D) 1 ), and by 
matching its asymptotic expansion at large z to the form it assumes for the spectral density 
of Eq. J7J). 



12 



ft(z) 




0123456789 1( 



z 

FIG. 1: Histogram: microscopic spectral density of the Dirac operator for a random matrix model 
with a color subblock and N c = 3. Solid line: spectral density predicted for xGUE (Nf — > 0). 
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FIG. 2: Histogram: microscopic spectral density of the Dirac operator for a random matrix model 
with a color subblock and N c = 2. Solid line: spectral density predicted for xGOE (Nj — > 0). 
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FIG. 3: Relative difference between the mean level spacing distribution for xGUE and that for 
a model with a color subblock and N c = 3. The errorbars are estimated from a series of nine 
independent runs of diagonalizations. 
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